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Abstract 

Recently we have shown that slow relaxations in the electron glass system can be understood 
in terms of the spectrum of a matrix describing the relaxation of the system close to a metastable 
state. The model focused on the electron glass system, but its generality was demonstrated on 
various other examples. Here, we study the noise spectrum in the same framework. We obtain 
a remarkable relation between the spectrum of relaxation rates A described by the distribution 
function P(X) ~ 1/A and the 1/f noise in the fluctuating occupancies of the localized electronic 
sites. This noise can be observed using local capacitance measurements. We confirm our analytic 
results using numerics, and also show how the Onsager symmetry is fulfilled in the system. 

PACS numbers: 71.23.Cq, 05.40.-a, 73.50.-h 
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I. INTRODUCTION 



1/f noise is ubiquitous and has been observed in a variety of physical systems, such 
as metals, semiconductors and superconductors 

HQ 

, as well as biological and economic 
systems p[ Q|. Systems displaying slow relaxations and aging are also abundant, and 



broad range of materia 
and biological systems 



Sj such as spin glasses jfi], structural glasses [6j, electron glasses 
have been shown to exhibit these properties. 
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There has been much interest in studying 1/f noise in electron glasses, both experimen 



tally and theoretically 
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22 1 has motivated 



the study of noise in electron glasses, particularly the study of the 'second noise' (describing 
the fluctuations in the noise spectrum), as a tool to distinguish between different glassy 
models. 

In a recent work 23|, we have proposed a model showing that the slow relaxations ob- 



served in the electrons glasses can be understood in terms of an underlying 1/ A distribution 
of relaxation rates, which arise naturally as the eigenvalues of a certain class of random 
matrices. It is the purpose of this work to show that the same model predicts a 1/f noise 
spectrum for the site occupancies, arising directly from the 1/A distribution. This estab- 
lishes a remarkable relation between 1/f noise and the universal slow relaxations. We study 
in detail the electron glass model, but the established connection between the two should 
be much more general than this specific system of interest, since the same equations can be 
used to describe a broad range of systems. 

The structure of the manuscript is as follows: we first define the model, and discuss briefly 
the local mean-field approximation used. We then review the previously derived results for 
the spectrum of the relaxation rates, and set up the equations describing the near-equilibrium 
fluctuations. We show that the Onsager symmetry is obeyed in the system, resulting in 
the specific form of the non-hermitian relaxation matrix, and show that using Onsager 's 
regression hypothesis we can relate the noise spectrum with the 1/A relaxation spectrum. 
Finally, we show numerical support for the calculation, and propose an experimental test 
for it. 
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II. DEFINITION OF THE MODEL AND THE LOCAL MEAN-FIELD APPROX- 
IMATION 



The model typically used to study electron glasses addresses a set of states, which are 
assumed to be Anderson localized due to the strong disorder present. While many studies 
have been performed, usually for numerical convenience, on a lattice, we choose to work with 
localized states positioned randomly in space, which we find more realistic. The number of 
electrons M is smaller than the number of sites N. The states' on-site energies contain disor- 
der, of typical magnitude W. Due to the localization, the Coulomb interactions between the 
electrons are not screened. They can not be neglected in the analysis and are not neglected 
within the model. The electrons are coupled to a phonon bath, assumed at thermal equilib- 



rium. Details of the model are found in 



24J , where we use a local mean-field approximation 



25| to study this model . We should emphasize that within this approximation the on-site 



and spatial disorder are not averaged over. We find that the approximation captures the 



2jj. In 



well-known results for the Coulomb gap [26|, |27| in the density-of-states (DOS), a soft gap 
emerging near the Fermi-energy as a result of the Coulomb interactions. It was also used 
to describe the slow relaxations experimentally observed in an aging experiment 
a later study we have shown that this model also describes the transition in the hopping 
conductance, from Mott's variable range hopping at higher temperatures to Efros-Shklovskii 
variable range hopping below the crossover 



An important point is that within the local mean-field approximation, the occupation 
numbers are continuous variables, and their dynamics is described by a coupled set of iV 
non-linear differential equations. Metastable states are configurations of occupation numbers 



for which the dynamics vanishes. In [24| we show that one can understand the dynamics 
near one of the many metastable states by linearizing the equations of motion close to it. 
This led us to an equation for the vector of deviations 5n from the equilibrium configuration 



in which the i'th site has the occupation n°: 

d5n 



dt 

where the explicit form of the non-diagonal elements of A is: 



A ■ Sn, (1) 
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with the equilibrium current from site z to site j (detailed balance is equivalent to 7 
being a symmetric matrix). From particle number conservation it follows that the sum of 
each column of A vanishes. Therefore the diagonal element of the matrix is minus the sum 
of the rest of the elements in the column, i.e., An = — Ylj^i Aij- 

Notice that exactly at the metastable state there would be absolutely no dynamics (or 
fluctuations) within this equation, even though it describes a system at finite temperature. 
This is clearly not physical, and in subsection IIV CI we add a necessary Langevin noise term, 
and determine its magnitude from thermodynamic considerations. 

III. KNOWN RESULTS FOR THE SPECTRUM 

It can be shown that the eigenvalues of the matrix A are real and negative. For a proof 
that the eig envalues are real see subsection IIVDI and for the negativity property see Ref. 



[17]. In 24j we have given numerical evidence that the model described in section HT1 yields 



a spectrum which is approximately: 

""Vrun 

where A max and A m i n are the upper and lower cutoffs of the distribution. In a later work 



291 ] we discuss a toy-model for this problem: we define a class of random matrices whose 
sum of columns vanish, and where the ij'th element decays exponentially with the distance 
of points i and j, distributed randomly. Notice that these are exactly the properties of the 
matrix A described in section [Til The merit of this simplification of the problem is that 
we manage to find an analytical formula for the averaged spectrum of the matrices, and 
verify the numerical result: up to logarithmic corrections in dimensions higher than one, the 
spectrum follows the 1/A spectrum in the low density regime. This is the essential result we 
would use in the following, and which would lead to the 1/f noise in the fluctuation of the 
occupation numbers. 
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IV. DERIVATION OF THE NOISE SPECTRUM 



Onsager's regression hypothesis, (see the Appendix for a derivation) states that the equa- 
tion of motion of the correlation function <f>ij(t) = (8ni(t)8rij(t)) is obtained by simply 
replacing <5n,(t) in the equation of motion ([I]) by the function 4>ij{t): 

^ = A ik ^{t). (4) 

To find the correlation function one also needs the initial conditions: 

<j> ij (0) = {6m(0)6n j (0))=l% 1 . (5) 

A. Equilibrium correlations 

By definition, the inverse of the matrix (3 describes the equilibrium correlations in the 
system. We would now like to calculate (3. Notice that generally /3 ij7 - = and that, for 
example, for a free Fermion system we have /3^j = n i (l — nj)Sij. 

In our case, we consider an interacting Fermion system, and it is not a-priori clear that 
we can use the above form. One might argue that within the local mean-field approximation 
the problem is essentially a non-interacting problem, with renormalized on-site energies due 
to the Coulomb interactions. It is therefore reasonable that the free Fermion result, not 
depending on the on-site energies, would still be valid. However, we shall now show that 
there is an additional term in j3 arising from the interactions. To do this, we write the free 



energy 



i i^j l 3 i 

with hi = Hi — ^. The local mean- field equations can be obtained from the minimalization 
condition = 0. Notice that each site contains a positive background charge of |, to keep 
charge neutrality. 

Expanding F near a metastable state (a local minima), we have: 

kT \ s 

F = F + — 22 PijSmSnj, (7) 

i,3 
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with: 



1 e 2 

^ = ^ n 0(i_ n 0) +^>-- ( 8 ) 
The correlation matrix is proportional to since we have a quadratic free energy: 



{5nM 3 ) = (9) 

The free fermion result is corrected for i ^ j by the interaction term. Notice that at low 
temperatures rij tend to or 1 exponentially, and thus the matrix is nearly diagonal, i.e., 
the interaction term is negligible. 

The matrix (3 has another physical meaning: if we introduce the conjugate variables to 
Hi, the so-called thermodynamic forces fiu defined as J-j^, we find that 

Hi = Pikn k . (10) 

B. Onsager symmetry 

According to Onsager 's principle (see Appendix) the kinetic coefficients, 

7« = ( n ) 

k 

are symmetric: 

Hj = la- (12) 
Indeed, Eqs. [2] and [8] are related according to: 

A = 7 /3, (13) 
where 7 is the symmetric matrix describing the transition rates at equilibrium, defined 



earlier. This puts the results of the linearization procedure of [24j in a much more natural 
context, and explains why the matrix A describing the relaxation is not hermitian in the 
general case: the multiplication of the symmetric matrix 7 and the symmetric matrix (3 is 
not expected to be (and is generally not) a symmetric matrix. 
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C. Langevin approach 



Due to the finite temperature, we expect the system to fluctuate around the equilibrium 
distribution. This can be modeled using a Langevin equation. In a similar fashion to the 
usual treatment of Brownian motion, where the equation governing the motion is % = 
— - + f(t), we add a Langevin noise term to Eq. ([!]): 

^ = A.5n + f(t), (14) 

where fi is a white-noise term. We shall now show, however, that the different components 
of the noise vector must be correlated, i.e., (fifj) = F is a non-diagonal matrix. As in the 
case of Brownian motion, thermodynamic considerations determine the magnitude of the 
noise, as well as the correlations between the different vector components. It is known that 



F = A/3~ l = 7 30f] . There is a clear physical intuition for this: correlations between sites i 
and j are contributed from the direct microscopic current between the two sites, 7^. Since 
7 is not diagonal, there exist correlations between the different components of the noise 
vector, as stated previously. See the Appendix for an explanation of the physical meaning 
of the matrix elements of 7 as kinetic coefficients. 

D. Form of the spectrum 

_____ ^2 

It is useful to diagonalize Eq. (TjJ by multiplying both sides on the left by ft'; . We then 



have, using ( 1131) and a matrix notation: 

&^ = ^^^ m . (15) 

Defining <fi = (5 1 l 2 <\) and 7 = /? 1 / 2 7/? 1//2 we reduce the equation to a symmetric form 
(notice that 7 = 7 T ): 

da) 

We can now diagonalize the symmetric operator 7 in a standard way: 

u^u = A (17) 

with A a diagonal matrix, and the elements of Uij = ip{ are the orthonormal eigenbasis 
(v)u = uv) = I) of 7. The eigenvalues of 7 turn out to be identical to those of A. Notice 



that this proves that the eigenvalues of A are real. To see this, it is useful to utilize the 
definition u~ l = v)f3 1/2 , u = (3~ 1/2 u =^> u = (3 1/2 u. Thus: 

A = M f /3 1/2 7/5 1/2 M = u-^^u = vT^Pu = u~ l Au (18) 

We can immediately write the solution for <fi: 

= e ~^4>(t = 0) = ue~ Kt u^(t = 0) => = /T^&r^/r 1 / 2 

Explicitly: 

^ = E A^ 1/2 C e " Aa W /2 - (19) 

This equation is an exact formula for the correlation function at site z, and will be 
later used to numerically find the spectrum. At low temperatures, as mentioned in section 
IIV A| the occupations tend to and 1 exponentially since they follow Fermi-Dirac statistics. 
For all sites whose distance from the Fermi-energy is larger than the temperature, we can 
therefore neglect the off-diagonal elements of the matrix. This is nevertheless a non-trivial 
approximation, since the sites close to the Fermi-energy are the ones contributing the most 
to the noise. For this reason we test the validity of this approximation numerically, and 
show that it doesn't change the form of the spectrum at the end of this section. Making the 
approximation we obtain: 

Pa 

a. 

Or in the frequency domain: 

This is an explicit formula for the noise spectrum. We can make a further approximation, 
by noticing that for most sites l//3u = rii(l — Hi) is very close to zero, since Hi is the 
Fermi function of the on-site energy, which is typically much greater than the temperature. 
Therefore, only a small fraction of sites for which the on-site energy (renormalized by the 
interaction) is of the order of kT away from the Fermi energy, will contribute to the sum. 
This is physically clear: it is exactly those sites which are close to the Fermi energy whose 
charge can fluctuate and contribute to the noise. Other sites are nearly permanently empty 
or full. For sites with energies Ei a distance much smaller than kT from the Fermi energy, 
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rij(l —rii) = 1/4. We therefore make an approximation and sum only over a partial number 
of the sites, for which we replace 1//^ by 1/4. This leads us to the equation: 



fcM~£>lf-^V < 21 > 

a a 

Taking another average over the different sites, we can replace (((VC) 2 )) by a constant 
jj, leading to: 

<<^~^£i+f^. (22) 

where ((; )) denotes averaging over sites as well as time. 

At this stage the equation takes the form leading to 1/f noise in various other theories: 
a discrete sum of equally contributing Lorentzians. The crucial thing is that the non-trivial 
weights arising from the matrix f3 did not affect the structure of the noise statistics. Finally, 



turning the sum into an integral and using the fact that P(\) ~ j in a large window 24j, 
we obtain: 

This shows that for \ m i n <C u> <C X max , a 1/f spectrum indeed follows for the noise in 
the average occupation number, which is one of our main results. In the future we intend 
to study the implications of the fluctuations in the occupation numbers on the conductance 
fluctuations. 

We have performed a numerical test of the calculation, taking the complete form of the 
matrix j3, including the off-diagonal elements, and using Eq. ( 1191) . Figured] shows the results 
for the spectrum, as well as explanations of the numerical procedure. The best linear fit 
over more than 4 decades gave a slope of —1.03 ± 0.03, close to the expected result. 



V. A PROPOSAL FOR AN EXPERIMENT 



20 



2l| 



A 1/f noise spectrum has been observed experimentally in electron glasses 
The measurements performed so far, however, were always for the conductance. This is not 

shows a strong connection between 



the quantity we calculate in this work, although 



the excess conductance and the deviations of the occupations. We hereby propose another 



10" 6 10" 5 10" 4 10" 3 10" 2 

CO 



FIG. 1: A numerical test of the noise spectrum. In the first step, the equilibrium occupation 
numbers, energies and transition rates were found for an arbitrary metastable state, using the 



procedure descrived in 24|. Next, the matrix (3 was constructed using Eq. (|8]). and the spectrum 
of A defined in Eq. (|13p was found. Finally, the spectrum was computed using the Fourier 
transform of Eq. (|19p . and averaged over the different sites. The spectrum was averaged over 100 
realizations. The best fit over more than 4 decades gave a slope of —1.03 ± 0.03. The system is 
one-dimensional, with iV* = 100 sites, and 20T = e 2 /r nn = W, where T is the temperature, r nn the 
average nearest-neighbor distance, and W the magnitude of the disorder. £/r n „ = 0.1, where £ is 
the localization length. 

possible experiment, in which the fluctuations in the occupation numbers are measured 
directly. In this experiment a small capacitor is placed in close proximity to the sample. The 
fluctuations in the occupation numbers of the neighboring sites would induce fluctuations in 
the voltage of the capacitor, which can be measured. Averaged over the different location in 
the sample, these should obey the 1/f noise discussed in section HVl A similar experimental 



technique has been implemented by 3l|, |32|, |33( . It would be very interesting to see whether 



a change in the nature of the noise is seen across the metal-insulator transition in such a 
measurement. 
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APPENDIX A: ONSAGER SYMMETRY OF KINETIC COEFFICIENTS 



It is the purpose of this Appendix to give a concise derivation of the Onsager regression 



30|. 



hypothesis and the Onsager symmetry of the kinetic coefficients, along the lines of 

Let us consider fluctuations of variables rij around some stable state in thermal equilib- 
rium. For convenience we take the equilibrium values to be 0. Close to the stable point 
we can expand the free energy to second order (the first order vanishes since we are at an 
extremum), and arrive at Eq. (|7j): 

F = 2^— PijUiUj. (Al) 
The Boltzmann distribution immediately gives us, after a gaussian integration, that: 

(^(0)71,(0)) =/^\ (A2) 

where () is the statistical ensemble average. We see that the inverse of (3 gives the 
equal-time correlation matrix. 

We now define the conjugate variable to n (the thermodynamic forces), as in section HVl 

dF 

dn = ^ ikUk ' 

A similar calculation to the one done for the correlations of n gives that (fiifij) = 
and (riiHj) = 5^. 

Let us assume the linearized equation of motion at the stable state is n = An. 
Defining (pij = (ni(t)rij(0)), we obtain that: 

fa = (riifijiO)) = A ik (n k (t) nj (0)). (A4) 
Thus we have the matrix equation: 

= A<j>. (A5) 

This is the Onsager regression principle, stating that the correlation function obeys the 
same relaxation equation as that of the microscopic variables. 
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Let us look at Cij(t) = = J" Assuming time-reversal symmetry , C is a 
symmetric matrix: Since we are at equilibrium we have time translation invariance, thus 
rfKOjn,, (o)) _ <iK(o)nj(-f)) ^ an( j us j n g the time-reversal symmetry this is equal to = 

Cji. Taking t = we obtain: 

C^(0) = (4^(0)71,(0)) = A iifc 3r//. (A6) 
Thus, A/5" 1 = 7 is a symmetric matrix. 

Notice that n = Af3~ l f3n = jjl, which is why the elements of 7 are called kinetic coeffi- 
cients, relating the fluxes (derivatives of n) to the forces \i. 
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